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ABSTRACT 

We present the barium surface abundance of 12 blue stragglers (BSs) and 18 main-sequence (MS) 
stars in the intermediate-age open cluster NGC 6819 (2.5 Gyr) based on spectra obtained from the 
Hydra Multi-object Spectrograph on the WIYN 3.5 m telescope. For the MS stars we find [Fe/H] = 
-1-0.05 ± 0.04 and [Ba/Fe] = —0.01 ± 0.10. The majority of the BS stars are consistent with these 
values. We identify five BSs with significant barium enhancement. These stars most likely formed 
through mass transfer from an asymptotic giant branch star that polluted the surface of the BS with 
the nucleosynthesis products generated during thermal pulsations. This conclusion aligns with the 
results from the substantial work done on the BSs in old open cluster NGC 188 that identifies mass 
transfer as the dominant mechanism for BS formation in that open cluster. However, four of the BSs 
with enhanced barium show no radial-velocity evidence for a companion. The one star that is in a 
binary is a double-lined system, meaning the companion is not a white dwarf and not the remnant of 
a prior AGB star. In this paper we attempt to develop a consistent scenario to explain the origin of 
these five BSs. 

Subject headings: open clusters and associations: individual (NGC 6819), stars: abundances 


1. INTRODUCTION 

Blue stragglers (BSs) are stars that do not follow the 
pathways laid out by standard stellar evolution. In star 
clusters they are typically identified as being brighter 
and/or bluer than the main-sequence (MS) turnoff. BSs 
have been identified in open and globular clusters, dwarf 
spheroidal galaxies, and the Galactic field. Gurrent mod¬ 
els invoke either mergers, collisions, or mass transfer 
events to explain their origins. 

Extensive work has been done on the BSs in the old 
open cluster NGG 188 (7 Gyr). Studies show that 80% of 
the BSs in NGG 188 are in binaries. Most have periods 
on the order of 1000 days and these long-period^ sys¬ 
tems have a secondary mass distribution peaked at 0. 5 
(iMathieu fc Gelleil I2009L iGeller fc MathieuI 1201111 . 
iGosnelll '( 201411 used the Hubble Space Telescope far- 
ultraviolet Advanced Gamera for Surveys/Solar Blind 
Ghannel to directly observe white dwarf companions of 
four BSs in NGG 188 and infer the presence of three 
more hot, young white dwarf companions (Teg ^ 11,000 
K, age < 400 Myr). All these lines of evidence indicate 
that mass transfer, particularly from an asymptotic gi¬ 
ant branch (AGB) star, is the dominant formation mech¬ 
anism for the BSs in NGC 188. 

Such mass transfer events can pollute the surface abun¬ 
dance of the BS with nucleosynthesis products from the 
evolved donor. Specifically, AGB mass-transfer should 
enhance the surface abundances of s-process elements, 
like barium, created during the thermally pulsing asymp¬ 
totic giant branch (TP-AGB) phase of stellar evolution. 
The other formation pathways, mergers and collisions, 
predict no such enhancements. This makes barium an 
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ideal marker for an AGB mass-transfer formation history 
of a BS. 

This same scenario (i.e., surface pollution from an 
AGB mass transfer event) is well understood to be the 
formation mechanism for barium stars. Barium stars 
are population I red giants that have enhanced abun¬ 
dances of s-process elements (particularly barium) and 
carbon, compared to other giants. They make up about 
1% of all G/K-type giants. On average they have sub¬ 
solar metallicities, masses between 1 Mm and 2.5 Me-, , 
and are in binaries (iJorissen et al.iri998LlPols et al.l[2003L 

lizzard. Dermine fc Churchll2010[l 

Previous efforts to use abundance anomalies to con¬ 
strain formatio n histories of BSs includes work on 
field BSs (e.g.. iSiyara rii et al.l 1200411 globular clusters 

/e.g..lFerraro et al.ll2006HLovisi et al.ll20i3aHLovisi et al.l 

2013bll. and opeii clusters (e .g., McGahee et al.l 2013 ; 
[Millir nan. Mathieu &: Schuled 1201311 . iSivarani et akl 
1^04ll analyzed a low-metallicitv halo BS in a spec¬ 
troscopic binary. They measured large overabundances 
of carbon, nitrogen, and oxygen as well as large en¬ 
hancements in s-process elements, and suggest that 
this star accreted m aterial from an AGB companion. 
iFerraro et al.l (1200611 discovered a subpopulation of six 
BSs in the globular cluster 47 Tuc that showed signif¬ 
icant depletion of carbon and oxygen com pared to the 
turnoff stars and the other BSs measured. iLovisi et al.1 
(|2013bll found a similar oxygen depletion in four BSs of 
M30. Both ascribe these depletions to the products of 
CNO burning being on the surface of the BS and con¬ 
clude these BSs to have a mass transfer origin with a red 
giant donor. 

For this paper, we focus on identifying an AGB mass 
transfer origin for BSs by measuring the barium abun¬ 
dance of BS and MS stars in the rich intermediate- 
age open cluster NGG 6819 (2.5 Gyr; a= 19^41™17®.5 
(J2000), 5= -b40°ll'47")- As one of four open clus¬ 
ters located within the original field of view of Ke- 
pler, NGG 6819 has experienced an increase in in- 
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terest in recent years. Fo r examp le, _ the WIYN 

Open Cluster Study"^ (WOCS; iMathied 1200011 h a s P^b - 
lished wide-held VI photometry (lYang et al.l I2013I1 . 
proper mo tions IPM. IPl atais et al.l 1201,111 . radial veloc¬ 
ities (RV, I Milliman et al.l 1201411 , and iron abundance 
measurements (fLee-Brown et al.M2015[l for this cluster. 
Other work has uncovered several peculiar objects in 
NGC 6819 i ncluding a candidate quiescent low-mass X- 
ray binary (iGosnell et al. 1201211 , possible evolv e d de¬ 
scendants of past BSs (jRosvick fc Vandenberd 119981 
iCorsaro et al.ll2 012ll. and a puzzling lithium-rich red gl - 
ant ( Anthonv-Twarog et al.ll20T^ iGarlberg et al.ll2015l l. 
The extensive information available for this cluster com¬ 
bined with its rich BS population made NGG 6819 an 
ideal target for our BS abundance study. 


2. TARGET SELECTION 

We selected targets for our abundance study that 
are identihed as m embers o f NGC 6819 based o n PM 
(IPlatais et al.ll2n i,‘lll and RV (I Milliman et'mi2014ll infor- 
mation. iMilliman et al.l (j201^ identihed 17 BSs in NGC 
6819 that are three-dimensional kinematic members, but 
that fall outside of the magnitude and color limits for 
normal stars and binaries from the cluster turnoff region 
or blue hook. We were able to obtain abundance mea¬ 
surements for 12 of these 17 BSs. Rapid rotation or hber 
separation requirements prevented us from getting mea¬ 
surements for the complete sample of BSs. To compare 
our BS sample to the standard abundance value in NGC 
6819 we also obtained measurements for 18 MS stars that 
are single cluster members with narrow spectral lines. 
Figure [I] shows the distribution of these targets in the 
color-magnitude dia gram (CMP ) for N GC 6819 using 
VI photometry from lYang et al.l (1201,811 . 

Table [T] lists the photometry, membership probabil¬ 
ities, e/i, and membership classihcation for the stars 
in our samp l e. Fo r a full explanation of eji please see 
iGeller et al'l (j2008ll . but briehy it is the ratio of the RV 
standard deviation (the external error, e) to the mea¬ 
surement precision (the internal error, i). The WOCS 
RV surveys use an eji threshold of 4 to distinguish be¬ 
tween velocity-variable stars and single members (SMs). 
The me mbership classification category is explained in 
detail in IMilliman et all (j2014fl and the distinctions are: 
cluster members have Prv ^ 50% and ^ 4%; SMs 
have eji < 4.0; binary likely members (BLMs) have an 
e/j > 4.0, but do not have a completed orbital solution; 
binary members (BMs) have an e/z ^ 4.0 and a complete 
orbital solution. IMilliman et al.l (1201411 find ^40% (7/17) 
of the BSs in NGC 6819 to be either BLMs or BMs. We 
are able to measure surface abundance values for five of 
these velocity-variable BSs. The two remaining binary 
BSs are too rapidly rotating (nsini = 60 km s“^ and 70 
km s“^) for us to get accurate abundance measurements. 

3. OBSERVATIONS AND DATA REDUCTION 

We obtained spectra of NGC 6819 in 2014 July 
with the Hydra Multi-object Spectrograph (jBarden et al.l 
1199411 on the WIYN 3.5 m telescope. Our observations 
included two different spectrograph setups covering two 
different wavelength regions centered at 5835 A and 6646 
A. The first setup has a wavelength range from 5700 to 


4 This is WIYN Open Cluster Study paper LXVIII 



Figure 1. CMD of all radial-velocity and proper-motion members 
of NGC 6819. To guide the eye a PARSEC zero-age main-sequence 
and a 2.5 Gyr isochrone are shown as the gray solid lines. Stars 
bluer than the dashed line are considered blue stragglers. Binary 
blue stragglers are marked by a circle and the one double-lined 
system is marked by two circles. Stars with measured abundance 
values are marked with large blue circles for the blue stragglers and 
large orange squares for our main-sequence control sample. 


5950 A with a resolution of i? ^ 19,000 and total ex¬ 
posure time of 5.5 hr. For these spectra we achieved 
a signal-to-noise ratio of 95 per resolution element for 
the faintest MS star {V = 15.29) and over 300 for the 
brightest BS {V = 1^82). The second spectral range 
extends from 6440 to 6840 A with a resolution of i? ~ 
13,000 and total exposure time of 8.7 hr. We achieved 
a signal-to-noise ratio in this spectral range of 150 per 
resolution element for the faintest MS star and over 415 
for the brightest BS. The results of this paper focus on 
the measured barium abundance from the Ba II feature 
at 5853.7 A and our iron abundance was based on Fe I 
lines from both spectral windows. 

The data were reduced using standard IRA F^ routines 
and fo llowing the basic procedure outlined in ISteinhaueil 
(|2003f) . In brief, the individual frames were trimmed 
and overscan subtracted before a master bias, made from 
combining a series of bias frames, was subtracted. Cos¬ 
mic rays were remove d from the object fr ames using the 
routine L.A. Cosmic (|van DokkumI[200lll . The spectra 
were then dispersion corrected, and the extracted spec¬ 
tra were flat-fielded, throughput corrected, and sky sub¬ 
tracted. Sky levels were determined from 27 sky fibers 
placed throughout the field of view that were averaged 
together to create a single sky spectrum. We corrected 
for Doppler shift using RVs derived from the centroid 

® IRAF is distributed by the National Optical Astronomy Ob¬ 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1 

Stellar Information and Parameters 


woes ID 


V 

V-I 

Prv 


eji 

Class 

'Te// (K) 

log ff 

Vt (km s 

1010 

BS 

12.82 

0.29 

86 

87 

8.66 

BM 

8524 

3.75 

1.0 - 4.0 

4004=^ 

BS 

13.98 

0.61 

94 

99 

17.45 

BM 

6900 

6650 

3.95 

4.0 

2.0 - 2.5 

2.0 

5006 

BS 

13.84 

0.58 

94 

99 

1.75 

SM 

6970 

3.75 

1.7 - 2.0 

8021 

BS 

14.18 

0.36 

93 

15 

1.29 

SM 

8006 

4.13 

2.0 - 3.0 

9003 

BS 

13.77 

0.69 

94 

99 

0.93 

SM 

6591 

3.60 

2.0 - 2.5 

10006 

BS 

14.37 

0.44 

79 

97 

1.98 

SM 

7419 

4.09 

2.0 - 3.0 

14015 

BS 

14.71 

0.50 

92 

97 

2.28 

SM 

7480 

4.14 

2.0 - 3.0 

16009 

BS 

14.79 

0.64 

93 

99 

5.72 

BLM 

6679 

3.97 

2.0 - 2.5 

16021 

BS 

14.72 

0.61 

84 

6 

1.94 

SM 

6820 

4.00 

2.0 - 2.5 

21022 


15.12 

0.77 

86 

47 

1.89 

SM 

6282 

3.96 

2.0 - 2.5 

22018 


15.03 

0.75 

86 

93 

0.43 

SM 

6352 

3.90 

2.0 - 2.5 

22022 


15.14 

0.74 

93 

53 

0.40 

SM 

6390 

3.96 

2.0 - 2.5 

23017 


15.11 

0.75 

94 

99 

1.19 

SM 

6353 

3.96 

2.0 - 2.5 

24027 


15.15 

0.79 

94 

58 

0.88 

SM 

6248 

3.96 

2.0 - 2.5 

25027 


15.26 

0.70 

93 

61 

1.09 

SM 

6508 

4.00 

2.0 - 2.5 

26017 


15.17 

0.72 

93 

99 

3.21 

SM 

6467 

3.96 

2.0 - 2.5 

26019 


15.21 

0.71 

94 

75 

2.83 

SM 

6549 

3.98 

2.0 - 2.5 

27019 


15.24 

0.73 

94 

98 

0.67 

SM 

6352 

4.00 

2.0 - 2.5 

27021 


15.18 

0.73 

84 

57 

1.24 

SM 

6389 

3.96 

2.0 - 2.5 

28017 


15.27 

0.71 

93 

99 

1.84 

SM 

6389 

4.00 

2.0 - 2.5 

30027 


15.29 

0.78 

57 

56 

0.32 

SM 

6317 

4.02 

2.0 - 2.5 

31028 


15.10 

0.75 

94 

60 

1.57 

SM 

6248 

3.96 

2.0 - 2.5 

32023 

BS 

15.28 

0.61 

94 

60 

4.73 

BLM 

6869 

4.19 

2.0 - 2.5 

32026 


15.11 

0.70 

93 

59 

0.73 

SM 

6428 

3.96 

2.0 - 2.5 

33038 


15.15 

0.78 

92 

36 

1.06 

SM 

6282 

3.96 

2.0 - 2.5 

41039 


15.11 

0.77 

94 

58 

0.37 

SM 

6352 

3.96 

2.0 - 2.5 

47030 

BS 

15.08 

0.62 

92 

21 

8.01 

BLM 

6820 

4.11 

2.0 - 2.5 

51059 

BS 

14.67 

0.66 

53 


1.19 

SM 

6635 

3.95 

2.0 - 2.5 

62034 


15.26 

0.71 

93 

47 

1.53 

SM 

6316 

4.00 

2.0 - 2.5 

63044 


15.29 

0.77 

85 

47 

0.27 

SM 

6123 

4.02 

2.0 - 2.5 

woes ID 4004 

is a double-lined system 

The 

spectral parameters for the 

primary 

are listed in 

the first entry and the secondary information is 


listed second. 


of a one-dimensional cross-correlation function with the 
fxcor task and an observed day-time sky spectrum. Indi¬ 
vidual exposures were added together to create the final 
spectra used in this analysis. 

4. ATMOSPHERIC PARAMETERS 

As is typical for abundance studies at the se spectral 
resolutions fe.g. lAnthonv-Twarog et ahll201Clll . we obtain 
atmospheric parameters based on photometric relation¬ 
ships and isochrones that we detail below. 

We determined the stellar parameters listed in Ta¬ 
ble [I]a|te£_ado£tin£_tlm reddening and distance modulus 
from lYang efaH (l20l3ll : E{B - V) = 0.14, E{V - I) = 
0.20, and {m — M)v = 12.36. Effective tempera tures 
were based on the BV color fr omiHole et al.l (1200911 and 
the fo llowing relationship from [Blackwell fc Lvnas-Gravl 

dUlH), 

Teff = [9134 - 8600(B -V) + 5398{B - Vf] 

' 0.40+ 1.30[Fe/H] l 

100 ■ ^ ^ 

We used the iBragaglia et al.l (l2001f l value of [Fe/H] = 
-1-0.0 9 for the metalli c ity. Adopting [Fe/H] = —0.02 
from ILee-Brown et al.l (120151) changes our temperatures 
by 1.2% or approximately 75 K for our MS stars. For 
the three stars without BV information (WOGS 33038, 
woes 41039, and WOGS 51059), temperature was 
based on stars nearby in the VI CMD that did have 
BV photometry. 

We adopted log g values from a 2.5 Gyr PAR- 


SEG (jBressan et al.l 1201211 isochrone for the MS stars 
and for the BSs we used PARSEC evolutionary tracks 
that passed through the L-Teff positions of these stars. 
We recognize this method for approximating BS grav- 
ity is simplistic based on uncertainties in BS masses 
([Mathieu fc Gelleill2015ll . but we take this into account 
by adjusting the log g values by ±0.2 dex in our analy¬ 
sis. We also note that [Fe/H] shows little dependence on 
log g and that microturbulence, Ut, has a much stronger 
effect on [Ba/Fe] than log g (see Table [2]). 

For the MS stars and fainter BSs we varied the micro¬ 
turbulence velocity from 2.0 km s~^ to 2.5 km s~^ based 
on the Vt results from lEdvardsson et alJ (|1993ll and 
iRamfrez. A llende P rieto fc LambertI ( 201311 . Based on 
the lTake^ et al.ri|2008| l study of A-tvne stars we chose a 
greater microturbulence range for the BSs in our sample 
with temperatures hotter than 7400 K. For WOGS 8021, 
woes 10006, and WOGS 14015 we varied Vt from 2.0 
km s“^ to 3.0 km s“^. For the hottest star in our sam¬ 
ple, woes 1010, we varied the microturbulence from 
1.0 km s“^ to 4.0 km s“^. BS WOGS 5006 has enough 
good quality iron lines that we could determine the range 
for Vt to be 1.7 km s“^ to 2.0 km s“^. The values for 
microturbulence that we used for each star are listed in 


Table [TJ 

woes 4004 is a special case in that it is a double-lined 
spectroscopic binary (SB2). The spectra and photometry 
have significant light contribution from both the primary 
and secondary star. For the atmospheric parameters we 
a dopted the primary an d secondary temperature values 
of iMilliman et al.] (|2014ll that were determined from cross 
correlating the spectra of WOGS 4004 against a grid of 
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synthetic spectra. Surface gravities were based on a com¬ 
parison between the temperature for each component, 
PARSEC evolutionary tracks, and the composite color 
and magnitude of WOCS 4004. We varied the microtur¬ 
bulence velocity from 2.0 km s“^ to 2.5 km s“^. 

In order to establish that our abundance results are 
significant and not highly dependent on our adopted pa¬ 
rameters we vary each quantity by an amount larger than 
the expected error. We shift Tes by ±200 K, log g by 
±0.2 dex, and Vt by ±0.5 km s“^ or more. The effects of 
each of these adjustments are discussed in more detail in 
Section [0] and listed in Tabled] 

5. ABUNDANCE ANALYSIS 

Our analysis was done using the 2014 J uly ver- 
sion of the LTE analysis code MOOG® (|Snedenl 
197311 incorporatin g the A TLAS9 model atmospheres of 
Castelli fc Kuruc j (1200411 . Our line list was compiled 
based on atomic information from the Vienna Atomic 
Line Database, (iKuoka et al.ll2000ll . We adjusted the 
log gf values slightly for two iron lines in order to better 
match our observed solar spectrum^ to expected abun¬ 
dance values. We also included hyperflne structure for 
the Ba II line at 5853.7 A as this better matched our 
day-time sky spectra with the expected solar abundance 
value. 

We measured the equivalent widths of our iron lines us¬ 
ing the spectrum analysis package SPECTRE®. We used 
these equivalent widths and the MOOG abfind routine 
to derive the abundance of each iron line. The [Fe/H]® 
abundance for each star was then determined by averag¬ 
ing the line-by-line differential abundances, i.e., the aver¬ 
age of all log V*(Fe l)iine i - log ^©(Fe ijune i for each 
star. The number of iron lines measured for each star is 
listed in Table |31 To determine the barium abundance 
we fit synthetic spectra to our observed spectra using 
the MOOG synth routine, focusing on the Ba II line at 
5853.7 A. Based on a solar model with atmospheric pa¬ 
rameters Teff = 5777 K, log g = 4.44, vt = 1.38 km s“^, 
and [Fe/H] = 0, we obtained log V 0 (Fe I) = 7.602 and 
log Al 0 (Ba) = 2.18. We determined the barium abun¬ 
dance relative to iron for each star using: [Ba/Fe]= (log 
VjBa) - logiy 0 (Ba|) - [Fe/H]*. 

iKorotin et al.l (1201 If) investigated NLTE effects on the 
[Ba/Fe] abundance of dwarf stars with similar temper¬ 
atures, gravities, and [Fe/H] to the stars in our sample 
and find a systematic increase of only ~0.05 dex. Based 
on this result we do not consider NLTE effects in this 
paper. 


5.1. The Double-lined System WOCS 4 OO 4 

WOGS 4004 is an SB2 with a period of 297 days, 
an eccent ricity of 0 . 21 , and a mass ratio, 7712 / 7711 , 
of 0.73 ([Miliiman et al'l 1201411 . We used a two- 
di mensional correlat i on tec hnique, TODCOR, developed 
bv iZucker fc Maz^ (jl994ll to determine the RVs of both 


® http://www.as.utexas.edu/~chris/moog.html 
^ For our solar spectrum we combined together day-time sky 
spectra taken with the same fiber configuration and reduced with 
the same methods as our science spectra. 

® http://www.as.utexas.edu/~chris/spectre.html 
® We use standard bracket notation to denote stellar abun¬ 
dances relative to solar values, e.g., [m/H] = log[A7(m)/7V(H)]* 
— log[A(m)/7V(H)]Q on a scale where log Af(H) = 12.0. 


the primary and secondary component. The two stars 
were separated by 30 km s“^ at the time of the ob¬ 
servation. We converted the spectra to a rest wave¬ 
length based on the RV of the primary star. As men¬ 
tioned in Section |4j the individual effective tempera¬ 
tures of the primary and secondary stars were derived 
from cross-correlation with a grid of synthetic spectra 
(iMilliman et al.l [201411 . The intensity ratio, ± 2 /^ 1 , of 
0.35 ± 0.04 was derived using a similar technique from 
spectra centered at 5187 A by Dr. G. Torres (2013, pri¬ 
vate communication). 

Because the absorption lines in the observed spectra 
have contributions from the primary and secondary stars 
determining [Fe/H] abundance from equivalent width 
measurements is problematic. Therefore we determined 
both [Fe/H] and [Ba/Fe] by fitting synthetic spectra to 
our observed spectra using the MOOG synth routine. 
For these fits we assumed the [Fe/H] values for the pri¬ 
mary and secondary to be equal. To obtain our abun¬ 
dance measurements we adjusted the iron and barium 
values of the synthetic spectra representing the primary 
and secondary stars until the combined spectrum best 
matched the observed spectrum as shown in Figured) 
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Figure 2. For the double-lined binary blue straggler, WOCS 4004, 
yye are able to use the pri mary and secondary temperatures derived 
bv IMilliman et ^ II2014I) . the luminosity ratio provided by Dr. G. 
Torres (2013, private communication), and TODCOR RVs to fit 
the observed spectrum (black dots) and derive the barium abun¬ 
dance for the individual stars. The best fit synthetic spectra (light 
blue) combines a hot primary (Teff = 6900 K, dark blue) with sig¬ 
nificant barium enhancement and a cooler secondary (Teff = 6650 
K, orange) with a barium abundance consistent with the expected 
cluster value. 



As described in Section [B] we repeated this fitting pro¬ 
cedure with a variety of temperature, gravity, and micro¬ 
turbulence values for the primary’s synthetic spectrum. 
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In order to limit the number of free parameters and be¬ 
cause the secondary star showed no unusual abundance 
features the atmospheric values for the secondary star 
were not adjusted. For each of the six best-fit synthetic 
spectra the iron and barium abundances of the secondary 
were near solar, but the barium abundance of the pri¬ 
mary was significantly above solar (see Table [3]) . 

6. ABUNDANCE RESULTS 

For our adopted atmospheric parameters, assuming a 
Vt = 2.0 km s“^, the MS stars in our sample have an 
average [Fe/H] = -1-0.05 ± 0.04 (the standard deviation 
of the distribut ion, s.d.). This value is consistent with 
the abundance iBraeaelia et al.l (|2001ll found from red 
clum p stars, [Fe/H] = -1-0.09 ± 0.03, and close to the 
value iLee-Brown et al.l (|2ni5h found from MS and turn¬ 
off stars, [Fe/H] = —0.02 ± 0.02. We find the average 
barium value for the MS sample is roughly solar, [Ba/Fe] 
= —0.01 ± 0.10 (s.d.). For the entire sample of BSs the 
average [Fe/H] = -1-0.06 ± 0.09 (s.d.). This is consistent 
with the above values, but with larger scatter perhaps 
caused by the higher temperature range of the BSs or 
because the atmospheric parameters are based on MS 
relationships and might be farther from the true BS val¬ 
ues. 

To ensure that our abundance results are robust and 
are not sensitive to our adopted atmospheric parameters 
we varied each parameter by an amount that exceeded 
the expected error. Specifically, for each star we shifted 
Teff by ±200 K, log g by ±0.2 dex, and varied the Vt val¬ 
ues as detailed in Sectional Table [3] outlines the average 
change in [Fe/H] and [Ba/Fe] from shifting the tempera¬ 
ture, gravity, and microturbulence values for the 18 MS 
stars in our sample. Temperature has the greatest im¬ 
pact on our [Fe/H] results, changing the abundance by 
a little over 0.1 dex. Our [Ba/Fe] results were most ef¬ 
fected by the change in microturbulence, with the 0.5 km 
s“^ increase resulting in 0.13 lower [Ba/Fe] values. 

With these changes in atmospheric parameters we mea¬ 
sured six [Fe/H] and [Ba/Fe] values for each star. The 
average of these six values is presented in Table IH along 
with the minimum and maximum [Fe/H] and [Ba/Fe] val¬ 
ues for each star. In Figure[3]we plot the average [Fe/H] 
versus average [Ba/Fe] for the MS and BS stars and use 
error bars to indicate the complete range in measured 
[Ba/Fe]. We also overplot the region corresponding to 
the 3a variation of the MS abundance values. 

As Figure [3] clearly shows, even allowing for consid¬ 
erable variations in temperature, gravity, and microtur¬ 
bulence, five BSs have barium abundances over the 3cr 
variation of the MS value for NGC 6819: WOCS 5006, 
woes 8021, WOCS 9003, WOCS 10006, and the pri¬ 
mary in WOCS 4004. These barium excesses are also 
evident in Figure [4] where we plot the observed spec¬ 
tra for four of these systems (see Figure [5] for the SB2, 
WOCS 4004). Along with our observed spectrum we 
plot two synthetic spectra with the atmospheric parame¬ 
ters of Table [TJ one has solar barium abundance and the 
other is our best fit spectrum that has significant barium 
enhancement. 

7. ABUNDANCE DISCUSSION 

These barium excesses indicate pollution of these five 
BS surfaces by an AGB companion and a mass transfer 



Figure 3. [Ba/Fe] vs. [Fe/H] for the stars in our sample. Blue 
stragglers are marked by blue circles and main-sequence stars are 
marked by orange squares. Binarity is marked by open circles and 
the primary of the SB2 WOCS 4004 is marked by two circles. The 
gray hatched box marks the 3a variation region of the MS values. 
The error bars mark the amount of variation in the [Ba/Fe] values 
from adjusting the temperature, log g, and vt values for all of the 
stars. 

origin for these systems, but there are some intriguing 
inconsistencies with this narrative that we explore below. 

7.1. Binarity of Barium-enhanced Systems 

The standard threshold used in WOCS RV surveys 
to determine if a star is velocity-variable is e/i > 4.0. 
However, e/i has a value less than 2.0 for all but one 
barium-enriched BS. With over 24 years of RV observa¬ 
tions these four low e/i targets are firmly categorized as 
non-velocity-variable. The one star that does have RV 
variation, WOCS 4004, has a companion that is defini¬ 
tively not a white dwarf and shows no means of providing 
the significant amount of barium measured on the sur¬ 
face of the BS. The current narrative of mass transfer 
from an evolved companion forming these BS systems, 
polluting the BS’s surface with barium and leaving be¬ 
hind a degenerate companion, does not evidently explain 
these systems. 

We suggest that the double-lined system, WOCS 4004, 
may have undergone a dynamical exchange after BS for¬ 
mation resulting in the current system with a luminous 
secondary. The normal abundance values of the sec¬ 
ondary is consistent with this scenario. It is extremely 
unlikely that all five of these enhanced systems under¬ 
went dynamical encounters after BS formation resulting 
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Table 2 

Abundance Shifts due to Atmospheric Parameters 



Uiff 

Tiff 

log g 

log g 



(-1-200 K) 

(-200 K) 

(-1-0.2 dex) 

(—0.2 dex) 

(-1- 0.5 km s“^) 

[Fe/H] 

-1-0.11 

-0.11 

-0.01 

-1-0.00 

-0.05 

[Ba/Fe] 

-0.02 

-1-0.04 

-1-0.08 

-0.07 

-0.13 


Table 3 

Abundance Information 

woes ID 


N 

[Fe/H] 

Range 

[Ba/Fe] 

Range 

1010 

BS 

3 

0.19 

-1-0.03 to -1-0.34 

-0.26 

-0.30 to -0.21 

4004^^ 

BS 


0.06 

— 0.05 to -1-0.10 

1.36 

4-1.10 to 4-1.50 

5006 

BS 

13 

0.09 

— 0.01 to 4-0.19 

0.94 

-HO.84 to -HI.16 

8021 

BS 

5 

0.20 

-|-0.08 to -|-0.33 

0.57 

4-0.42 to 4-0.67 

9003 

BS 

7 

-0.06 

-0.16 to -h0.04 

0.48 

4-0.27 to 4-0.59 

10006 

BS 

7 

0.16 

-|-0.0T to -|-0.2T 

0.99 

4-0.62 to 4-1.14 

14015 

BS 

3 

0.03 

— 0.06 to -1-0.13 

-0.14 

-0.22 to -0.07 

16009 

BS 

6 

-0.00 

— 0.10 to 4-0.10 

-0.11 

—0.18 to —0.03 

16021 

BS 

3 

-0.01 

-0.11 to -hO.lO 

0.24 

4-0.10 to 4-0.35 

21022 


4 

0.06 

— 0.05 to 4-0.19 

-0.12 

-0.22 to -0.02 

22018 


5 

0.02 

— 0.08 to 4-0.14 

-0.18 

—0.29 to —0.09 

22022 


8 

0.06 

— 0.04 to 4-0.18 

0.06 

—0.06 to -HO. 17 

23017 


7 

0.07 

-0.03 to -h0.20 

-0.11 

-0.21 to -0.03 

24027 


5 

0.05 

— 0.06 to 4-0.17 

-0.11 

-0.20 to -0.04 

25027 


5 

0.02 

— 0.08 to 4-0.13 

0.15 

4-0.01 to 4-0.24 

26017 


3 

0.07 

-0.04 to -h0.19 

0.02 

—0.09 to 4-0.12 

26019 


5 

0.01 

— 0.09 to 4-0.12 

-0.08 

-0.17 to -0.02 

27019 


8 

0.04 

— 0.06 to 4-0.16 

-0.05 

—0.14 to -HO.05 

27021 


7 

0.00 

-0.10 to -h0.12 

0.01 

—0.11 to 4-0.09 

28017 


8 

0.07 

— 0.03 to 4-0.18 

0.03 

-0.10 to +0.13 

30027 


6 

0.03 

— 0.08 to 4-0.15 

0.06 

—0.06 to +0.17 

31028 


11 

-0.02 

-0.12 to -hO.lO 

0.07 

—0.04 to +0.18 

32023 

BS 

5 

0.05 

-0.04 to -h0.15 

-0.00 

—0.07 to +0.07 

32026 


6 

0.05 

— 0.05 to 4-0.17 

-0.03 

—0.13 to +0.06 

33038 


11 

0.09 

-0.02 to -h0.20 

-0.16 

—0.26 to —0.08 

41039 


5 

0.11 

4-0.01 to 4-0.23 

0.09 

-0.04 to +0.20 

47030 

BS 

4 

0.00 

-0.09 to -HO.11 

-0.03 

-0.10 to +0.04 

51059 

BS 

6 

-0.04 

-0.14 to +0.07 

-0.01 

—0.08 to +0.08 

62034 


3 

0.01 

-0.10 to -HO.14 

0.07 

—0.06 to +0.18 

63044 


8 

-0.06 

— 0.16 to 4-0.07 

-0.14 

—0.22 to —0.05 


® woes ID 4004 is a double-lined system. The abundance values for the primary are listed. 


in four non-velocity variable stars and one double-lined 
binary. Especially since dynamical interactions for bi¬ 
naries with periods under 10^ days typical ly result in 


_ itYE _ 

shorter-period, tighter orbits (|Heggiel Il975ll . With no 
other mechanism available to produce the barium abun¬ 
dances observed besides an AGB companion and the low 
probability of the dynamical destruction of all these bi¬ 
naries, we look to other means of explaining these results, 
none of which offers a completely satisfactory solution. 

One possibility is that the four enhanced systems with 
no RV variation have companions that we are not able 
to detect. We modeled this scenario with the same 
method used to determine the overall binary detec- 
tion completeness of the WOCS RV survey outlined in 
Geller fc Mathieul (1201211 and detailed for NGC 6819 in 


Milliman et al.l ( 2014f) . Specifically, we generated ar¬ 


tificial binaries with period and eccentricity distribu- 
tions of the Ga l actic field binary population found by 
iRaghavan et al.1 (l201fll l. We assumed a secondary white 
dwarf mass of 0.6 Mq and a primary mass that ranged 
from 1.6 Mq to 2.4 Mq. Orbital inclination, longitude 
of periastron, and phase were chosen randomly. We used 
the actual observation dates and measurement precisions 
in the model. RV measurements were produced for a 


large number of synthetic binaries and their detectabil¬ 
ity assessed. Our detection completeness depends most 
strongly on the period of the binary system and has lit¬ 
tle dependence on the mass ratio. For each mass ratio 
our detection capability falls below 50% for periods over 
15,500 days and below 25% for periods over 40,000 days. 

Such long periods severely restrict the amount of mass 
that can be transferred to a companion. Adjusting pa- 
ram eters within the Binary-stell ar Evolution code (BSE) 
from iHurlev. Tout &: Pols((l2002D . the maximum mass we 
are able to transfer that resulted in a binary orbit over 
15,500 days was under 0.12 Mq. This value required 
highly eccentric systems (e ^ 0.5). 

A similar low mass e nhancement was also suggested 
by iGosnell et al.1 (1201411 for a long-period (3030 ± 70 
days) BS with a detected white dwarf companion in NGG 
188. Using BSE, they found the accretor received 0.101 
Mq through wind accretion. Surprisingly, this BS is one 
of the m ost luminous BS bina ries in NGG 188. As men¬ 
tioned in IGosnell et al.l (|2014ll this raises many questions 
about the mass-luminosity relation in BS stars, whether 
wind accretion is the correct scenario, and if the AGB 
wind prescription used in BSE is correct. The forma¬ 
tion of these four potentially very long period barium 
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Figure 4. The spectra for the blue stragglers with significant barium excess are shown (black dots). Overplotted are the computed profiles 
with [Ba/Fe] = 0.0 (red line) and the fitted barium abundance (blue line). 


enhanced stars adds to these questions. 

Another explanation for these non-velocity variable 
stars arises from work done with barium stars, red gi¬ 
ants with enhanced carbon and s-process ele ments. Bi¬ 
nary popu lation synthesis studies including iPols et ahl 
lj200.‘ll ) and lizzard. Dermine fc Churchl (|2010f ) seek to ex¬ 
plain the origin of barium stars and match the observed 
population fraction, measured overabundances, and the 
distribution of periods and eccentricities of these mostly 
binary systems. Without modified physics these stud¬ 
ies do not reproduce the observed barium star period- 
eccentricity distribution. Specifically they produce no 
eccentric systems with P < 3000 days and generate too 
many long p eriod (P > 10^) barium stars. In order to 
address this, lizzard. Dermine fc Churchl (|2010D apply a 
white-dwarf kick in their model as a way to generate the 
eccentricities observed in binary barium stars (although 
the mechanism for producing the white dwarf kick is un¬ 
known) . 

In the model this white dwarf kick also pushes weakly 
bound, long-period systems into wider orbits and even 


disrupts some binary systems. Either of these two sce¬ 
narios could explain the non-velocity-variable enhanced 
systems seen in NGC 6819. Unfortunately observations 
of actual barium stars rule out a large population of long 
period systems and the overwhelmingly large fraction of 
barium stars detected in binarie s means a very small dis¬ 
ruption rate for these systems (j.Torissen et al.l[T9^ . so 
a strong white dwarf kick seems unlikely. 

It is unclear why, if these enriched open cluster BSs and 
field barium stars all formed through the same process, 
their orbital properties would be so different. Hopefully 
more complete binary evolution and A^-body modeling 
for NGC 6819 will provide some insight. 

7.2. Amount of Barium Enhancement 

The five barium-enhanced BSs have [Ba/Fe] ranging 
from -kO.48 (WOCS 9003) to -kl.36 (WOCS 4004). We 
have suggested that the origin of the excess barium on 
these BSs is the s-process element production that hap¬ 
pened in the interior of AGB stars and subsequently was 
mass transferred to the BSs. Here we address whether 
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these levels of barium enhancement are feasible. 

First we assume that the AGB companion evolved as 
it would in isolati on. This is appropriate because, as 
discussed in section mi if these non-velocity-variable BS 
with barium excess are in binaries they have long-periods 
and the binary would not have significantly disturbed the 
stellar evolution, thermal pulses, or barium production 
of any AGB companion. 

We also assume that an accreting star would have a 
mass near the turnoff mass of NGG 6819, ~1.5 Mq, 
and would not have efficient surface convection zones 
that would significantly dilute the barium abundance 
of any accreted mat erial. This is supported by MESA 
(jPaxton et al.l 120111) models of solar metallicity stars, 
where a 1.4 Mq and 1.2 Mq star have only ^0.0006 
Mq and 0.0015 Mq^ respectively, in a surface convection 
zone at an age of 2.5 Gyr. 

The accreted material would have a higher mean 
molecular weight than the surface layer of the ac- 
cretor and thermohaline mixing should eliminate the 
molecular weight difference on thermal t ime-scales 
(|KiDDenhahn. Ruschenulatt fc ThomasI Il98flh . How¬ 
ever, abundance measurements of carbon-enhanced 
metal-poor stars have called into question the ef¬ 
ficiency of thermohaline mixing and vari ous reduc¬ 
tion scenarios have been investigated (lAoki et ah] 
20081 IStancliffe &: Glebbeekll2008l : [Thompson et al.ll2008t 

Bisterzo et al.l[2nilj ). Based on this we assume no signifi¬ 

cant amount of dilution of the barium surface abundance 
has taken place since the time the enriched material was 
accreted. 

T o compare our observ ational results with theory we 
use iGristallo et al.l (l2011[) models of the nucleosynthesis 
processes of AGB stars of different masses. These results 
are available from the FRUITY Database^° that provides 
the isotopic composition from H to Bi as it changes af¬ 
ter each third dredge-up (TDU) episode and the stellar 
yields produced for each model. Kepler asteroseismology 
puts the ave rage mass of red c lump stars in NGG 6819 
at 1.64 Mq (jMiglio et al.ll20l'^ . The FRUITY database 
has models for 1.5 Mq and 2.0 Mq. Assuming solar 
metallicity (z = 0.014) and no initial rotational velocity, 
the 1.5 Mq model experiences 5 TDU episodes with a 
final [Ba/Fe] = -1-0.65 and the 2.0 Mq model undergoes 
9 TDU events with a final [Ba/Fe] = -1-1.03. 

All of this indicates that the amount of barium pro¬ 
duced by a ~1.5 Mq to ^2 Mq initial-mass AGB com¬ 
panion can broadly account for the amount of barium 
excess seen in these NGG 6819 BSs. 

7.3. Age and Position in the CMD 

As outlined in the previous section, the levels of barium 
excess for the three stars with the highest enrichment 
(WOGS 4004, WOGS 5006, and WOCS 10006) require 
a donor with ~2 Mq initial mass. Based on the 1.5 
Gyr evolutionary lifeti me of a 2.0 Mq star (from the MS 
through the TP-AGB; iBressan et al.i [21)1^ and the 2.5 
Gyr cluster age of NGG 6819 that places the formation 
of these BSs at roughly 1 Gyr ago. 

The barium enhancement levels of the other excess BSs 
(WOGS 8021 and WOCS 9003) are more in line with the 
expected [Ba/Fe] levels of the current AGB stars at 1.64 

http://fruity.oa-teramo.inaf.it/ 


Mq ([Miglio et al.ll201^ . So we expect these BSs to be 
young systems and to have ages much less than 1 Gyr, 
but we are unable to quantify when they formed. 

We can compare the age esti mates of ■^l Gyr to the so¬ 
phisti cated A^-body models of iGeller. Hurley fc Mathi^ 
(1201311 for NGG 188. These models track BS formation 
throughout the cluster’s lifetime up to the present NGG 
188 age of 7 Gyr. Taking snapshots of the 20 realiza¬ 
tions of the NGG 188 model at a cluster age between 
2 and 3 Gyr, we construct the age distribution of mass 
transfer-formed BSs (i.e., the time since the end of mass 
transfer) that we expect in a cluster like NGG 6819. We 
plot this age distribution in Figure [S] and show that at 
a cluster age of 2—3 Gyr the mass transfer-formed BSs 
are typically young systems, but there is a long tail that 
extends to roughly 1.6 Gyr. Our Monte Carlo-sampling 
of this distribution shows that finding three or more BSs 
with ages older than 1 Gyr is not restrictively uncom¬ 
mon, happening ^20% of the time when we assume a 
total of 12 mass transfer-formed BSs. 

These three BS age estimates of ~1 Gyr are much older 
than the seven other known open cluster BS ages that 
range from 70 Myr to 400 Myr. These other BS ages 
were all determined fo r BSs i n the old open c luster NGG 
188 bv iGosn^ (|2014ll . The iGosnelll (120141) study was 
designed to find the most recently formed mass transfer 
BSs through detection of the UV excess produced by hot, 
young white dwarf companions. On the other hand our 
study determines age estimates for only the most barium 
rich and oldest BSs. So we consider these different age 
determinations to be a consequence of the observational 
biases of the two surveys. 



Figure 5. Histogram of the ages for mass transfer-formed BSs 
from the NGC 188 T^-body model of IGeller. Hurley fc MathieuI 
ll20il) from when the cluster was 2—3 Gyr. 


We highlight the CMD position of the five NGC 6819 
BSs with significant barium excess in Figure [51 All five 
BSs occupy a narrow region and fall between the evo- 
lutionary tracks of a ^1.75 Mq and a '^1.95 Mq star 
([Bressan et al.l[20T^ . We do not make mass estimates for 
these stars, because it is not clear if th ese standard tracks 
are ap propriate for BS mass estimates (jMathieu fc Gell^ 
120151) . but this positioning would suggest that these 
BSs have similar masses. The barium excess systems 
have varying distances from the zero-age main-sequence 
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(ZAMS) which might be attributable to age differences. 
However, their positions have no apparent dependence 
on barium abundance values. The two BS closest to the 
ZAMS are WOCS 8021, [Ba/Fe] = +0.57 and WOCS 
10006, [Ba/Fe] = +0.99. The three furthest from the 
ZAMS are the primary of WOCS 4004^^, [Ba/Fe] = 
+ 1.36; WOCS 5006, [Ba/Fe] = +0.94; and finally WOCS 
9003, [Ba/Fe] = +0.48. This distribution shows that the 
most barium-enhanced, and presumably oldest BSs at 
'^l Cyr, are not the farthest from the ZAMS. 

If in fact all of these BSs have comparable masses, 
these results suggest that barium abundance or CMD 
position are not accurate measures for BS age. Addi¬ 
tional evidence that BS age and CMD position is not 
a straight-forw ard relationship comes from NCC 188. 
iCosnelll (|2014ll found that the distance from the ZAMS 
is not necessarily equivalent with BS age. This conclu¬ 
sion was based on the positions off the ZAMS of the 
definitively very young BSs in NCC 188. iCosne n il2nTii 
suggests that BS ages based on single-star isochrones are 
unreliable. Our results also seem to suggest this although 
this conclusion is clouded by our assumptions. 

As for why the barium-enhanced systems in NCC 6819 
occupy such a narrow region in the optical CMD, we have 
no satisfying explanation. 



Figure 6. Same as Fig. [T] except BSs with significant barium 
excess are marked by red diamonds. 


8. FORMATION MECHANISMS FOR THE OTHER BLUE 
STRAGGLERS 

This order is based on the expected position of the primary 
from a simplistic deconvolution of the SB2 photometry. 


The origin of the other 7 NCC 6819 BSs in this abun¬ 
dance survey is still an open question. The lack of signifi¬ 
cant barium enhancement itself does not rule out a mass- 
transfer formation history for a BS. If the BS progenitor 
binary had a close enough orbit, the mass transfer event 
could have occurred before the giant companion made 
it to the ACB phase, or the binary could have stripped 
enough of the ACB envelope to prevent thermal pulses 
and any substantial barium production. Future work will 
focus on other abundance tracers that are indicative of 
earlier mass transfer activity. We will particularly fo¬ 
cus on oxygen and carbon anomalies that would indicate 
that mass transfer originated along the re d giant branch 
as seen in some g]obular cluster BSs (e.g.. iFerraro et aP 
120061 : iLovisi et aHl2013bll . 

It is harder to explain the lack of barium enrichment 
for long-period BS binaries. WOCS 1010 has an orbital 
period of 1,144 days, an eccentricity of 0.55, and has the 
lowest [Ba/Fe] of any of the BS surveyed. WOCS 1010 is 
also the hottest BS in NCC 6819 with a temperature of 
8524 K, meaning its surface abundance is most likely 
strongly affected by ra diative levitation (|Lovisi et al.l 
l2012LlLovisi et al.]l2013atl and we are unable to draw any 
solid conclusions from the abundance results. 

Another long period BS binary is WOCS 47030. It 
has a Teff = 6820 K and does not have a completed or¬ 
bital solution, but the RV data suggest an eccentric or¬ 
bit (e ^ 0.7) with a period around 3700 days. At that 
period the most likely mass transfer formation scenario 
would involve wind accretion from an ACB companion, 
but no enrichment in [Ba/Fe] is measured. This leaves 
open the possibility that WOCS 47030 formed from the 
merger of a hierarchical triple system as described in 
iPerets fc Fabryck"^ (|2009l l _ 

We note th at although iCeller fc MathieiJ (|20lHl and 
ICosnelll (|2014f l establish mass transfer to be the dominant 
BS formation process in NCC 188 and we find at least 
^30% (5/17) of the BS in NCC 6819 to be formed from 
specifically ACB mass transfer events, we still expect 
collisions and mergers to have formed some BSs in NCC 
6819. 

9. SUMMARY AND DISCUSSION 

As a way to distinguish between various BS formation 
mechanisms we measured the abundance of the s-process 
element barium in 12 BSs in the intermediate-age open 
cluster NCC 6819. Our observations of these BSs and 18 
MS stars were made in 2014 July using the Hydra Multi¬ 
object Spectrograph on the WIYN 3.5 m telescope. 

For our control sample of MS stars, we find an average 
[Fe/H] = +0.05 ± 0.04 and [Ba/Fe] = -0.01 ± 0.10. This 
iron abundance is consistent with previous spectroscopic 
abundance studies of NCC 6819. The average [Fe/H] 
value for the BSs is also consistent, [Fe/H] = +0.06 ± 
0.09. 

Out of the 12 BSs in our sample we identify five that 
have significant barium enrichment. From the amount of 
barium enhancement on these five BSs we find that: 

1. A ^1.5 Mq to ~2.0 Mq initial-mass ACB com¬ 
panion is able to produce that amount of [Ba/Fe] 
measured. 

2. This ACB mass range matches well with the cur¬ 
rent NCC 6819 red clump mass determined from 
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Kepler observations. 

3. The three most enriched BSs formed roughly 1 Gyr 
ago. 

We suggest that these barium enhanced systems 
formed via mass transfer from an AGB binary companion 
that polluted the BS with the nucleosynthesis products 
of the AGB star. This formation scenario is complicated 
by the fact that: 

1. Four of the enhanced systems show no RV evidence 
for a companion. 

2. If these non-velocity-variable systems are indeed in 
long-period binaries, we place the expected orbital 
periods at over 15,500 days. 

3. Such wide orbits severely limits the amount of mass 
that can be transferred to a companion and we find 
the maximum amount transferred to be under 0.12 

Mq. 

4. The one enriched system in a binary has a luminous 
companion with no means of providing the excess 
barium observed. 

5. This proposed formation scenario is similar to how 
barium stars form, but we find the barium-enriched 
BSs in NGG 6819 have very different binary prop¬ 
erties than these enhanced field giants. 

Future work should focus on incorporating more de¬ 
tailed binary evolution modeling for these barium en¬ 
riched BSs, perhaps with a different AGB wind prescrip¬ 
tion, and A^-body models specifically for NGG 6819 to 
help explain the origin of these enriched systems and 
their lack of RV-detected companions. Additional abun¬ 
dance tracers, such as other neutron capture elements 
and carbon, may also shed light on the formation histo¬ 
ries of these enhanced BSs, as well as the other BSs in 
NGG 6819. 

Moving beyond NGG 6819, we predict that there is 
an upper limit on cluster age for finding barium en¬ 
hanced BSs formed via mass transfer from an AGB bi¬ 
nary companion. The FRUITY database has no TDU 
episodes or s-process element production occurring for 
a solar-metallicity 1.3 Mq AGB star. Presumably this 
is because mass loss has reduced the envelope mass be¬ 
low the ~0.4—0.5 thr eshold required for TDU to 

occur (|Cristallo et aLllMIlll . We estimate this barium- 
enhancement cluster age limit by combining the standard 
stellar lifetime of a 1.5 Mq star (the lowest mass model 
on the FRUITY database that produces barium) and 
the amount of time a star “survives” as a BS. We use 
a PARSEC evolutionary track to get a standard stellar 
lifetime of 3 Gyr. The tail of the Wbody mass transfer- 
forme d BS age distribution of IGeller. Hurley fc Mathi^ 
(|2013ll provides us with ^2 Gyr as the typical upper end 
of BS survival time. So we place this cluster age limit 
at roughly 5 Gyr and do not expect barium-enhanced 
BSs in open clusters^^ older than 5 Gyr. Indeed, the 
surface abundance study of the BSs in the 7 Gyr open 

Lower metallicity environments will have a different age limit. 


cluster NGG 188 found no obvious b arium enhancements 
(|Milliman. Mathieu fc Schulerl[2M^ . 

We recommend that future BS abundance surveys look 
for s-process element enhancements in intermediate- 
aged open clusters younger than NGG 6819, but we do 
caution that clusters that are too young will potentially 
have more high temperature BSs affected by radiative 
levitation. 
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